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ABSTRACT 

We present moderate-resolution spectroscopic observations of red supergiants (RSGs) in the low- 
metallicity Local Group galaxies NGC 6822 (Z = 0.4Zq) and Wolf-Lundmark-Melotte (WLM; 
Z = O.IZq). By combining these observations with reduction techniques for multislit data reduc- 
tion and flux calibration, we are able to analyze spectroscopic data of 16 RSGs in NGC 6822 and 
spectrophotometric data of 11 RSGs in WLM. Using these observations we determine spectral types for 
these massive stars, comparing them to Milky Way and Magellanic Clouds RSGs and thus extending 
observational evidence of the abundance-dependent shift of RSG spectral types to lower metallicities. 
In addition, we have uncovered two RSGs with unusually late spectral types (J000158. 14-152332.2 
in WLM, with a spectral type of MS I, and J194453.46-144552.6 in NGC 6822, with a spectral type 
of M4.5 I) and a third RSG (J194449.96-144333.5 in NGC 6822) whose spectral type has varied 
from a M2.5 in 1997 to a K5 in 2008. All three of these stars could potentially be members of a 
recently-discovered class of extreme RSG variables. 



1. INTRODUCTION 

Red supergiants (RSGs) are the evolved He-burning 
descendants of moderately massive (1O-25M0) OB stars. 
Although RSGs are not the most massive or most lumi- 
nous members of the massive star population, they are 
the coolest and largest (in physical size), and as such rep- 
resent an important extreme in the H-R diagram (HRD) 
that has, in the past, been poorly matched by evolution- 
ary tracks. 

The large radii and cool temperatures of RSGs are con- 
strained by the Hayashi limit, the coolest temperature 
at which a fully convective star would remain hydrostat- 
ically stable (Hayashi & Hoshi 1961). The location of 
the Hayashi limit is dependent on metallicity, shifting 
to warmer temperatures (and therefore earlier spectral 
types) at lower metallicities (Elias et al. 1985, Levesque 
et al. 2006). However, for many years the placement of 
the Hayashi limit on the HRD was at odds with the "ob- 
served" locations of RSGs on the HRD, with the stars 
appearing to be considerably cooler and more luminous 
than theoretically permitted (Massey 2003, Massey & 
Olsen 2003). 

In recent years, we have used the new generation of 
MARCS stellar atmospheres to fit moderate-resolution 
spectrophotometry of RSGs in the Milky Way (Levesque 
et al. 2005), the Magellanic Clouds (Levesque et al. 
2006), and M31 (Massey et al. 2009). This work has 
resulted in a significant revision of these stars' effective 
temperature scales and their position on the H-R dia- 
gram. For the Milky Way and the Large Magellanic 
Cloud, the location of the RSGs is now in excellent agree- 
ment with the predictions of the evolutionary tracks. We 
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also see that the distribution of spectral types in these 
galaxies shifts with metallicity, yielding earlier average 
RSG spectral types in lower-metallicity galaxies. While 
this shift is partially due to lower Ti abundances yield- 
ing weaker TiO band strengths, it can primarily be at- 
tributed to RSG populations in lower-metallicity galaxies 
having warmer median effective temperatures. This is a 
consequence of the metallicity-dependent Hayashi limit, 
in accordance with the predictions of evolutionary mod- 
els (e.g. Schaerer et al. 1993, Charbonnel et al. 1993, 
Maeder & Meynet 2001, Meynet & Maeder 2005). The 
same effect is also observed in the red giant branches of 
globular clusters: a decrease in the number of heavy ele- 
ments contributing electrons corresponds to a decrease 
in opacity and a corresponding increase in the stars' 
surface temperature (Elias et al. 1985, Levesque et al. 
2006). However, the lowest metallicity RSGs in our sam- 
ple, from the SMC, remain slightly cooler than the evo- 
lutionary models predict, with some RSGs in the sample 
remaining to the right of the Hayashi track. It is unclear 
what exactly is causing this discrepancy. What is clear 
is that massive stellar evolution at low metallicity, par- 
ticularly in the later stages of these stars' lives, is still 
poorly understood. 

Additionally, three RSGs in the SMC and one in the 
LMC were found to display unusual variable behavior 
that has not been previously associated with RSGs, all 
showing considerable variations in V magnitudes and ef- 
fective temperatures (and spectral types) (Levesque et al. 
2007, Massey et al. 2007a). These stars suffer dramatic 
physical changes on timescales of months - when they 
are at their warmest they are also brighter, more lumi- 
nous, and show an increased amount of extinction. The 
extinction is characteristic of the effects of circumstellar 
dust (Massey et al. 2005) and is thought to indicate spo- 
radic dust production from these stars in their coolest 
states. One parameter these variable stars share is un- 
usually late spectral types; two of these stars, HV 11423 
and [M2002] SMC 055188, were previously observed in 
an M4.5 I state, considerably later and cooler than any 
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other star in the SMC and lying well to the right of the 
Hayashi limit. The physical explanation for these stars' 
rapid variations and unusual parameters, and their con- 
nection with low-metallicity environments, still remains 
to be explained. 

The metallicities of the RSG populations studied in 
Levesque ct al. (2005, 2006) spanned a factor of 5. By 
comparison, galaxies in the Local Group with known 
RSG populations span a factor of 15 in metallicity 
(Massey et al. 2007b). M31 represents the high metallic- 
ity extreme, with some studies placing its metallicity as 
high as 2x solar (Blair et al. 1982, Zaritsky et al. 1994; 
however, sec Massey et al. 2009 for further discussion 
of the M31 abundance). At the other extreme lies the 
Wolf-Lundmark-Melotte (WLM) galaxy, which, at 0.1 x 
solar (Hodge & Miller 1995) has the lowest metallicity 
of any Local Group galaxy currently forming stars (van 
den Bergh 2000, Massey 2003). Bresolin et al. (2006) 
have previously obtained spectra of 34 carly-typc super- 
giants in WLM, but spectra of evolved supcrgiants in this 
galaxy have never before been observed (nicinbcrship of 
the four suspected G supergiants detected in Bresolin et 
al. 2006 is not known for certain; for further discussion 
see Massey et al. 2007b). 

Here we present moderate-resolution spectroscopy of 
16 RSGs in NGC 6822 and 11 RSGs in WLM. The stars 
were selected using UBVRI from Massey et al. (2007b) 
and observed using IMACS on the Magellan 6.5m Baade 
telescope (§ 2). We determined spectral types for the 
sample (§ 3), comparing these stars to the Milky Way 
and Magellanic Cloud RSGs and extending observations 
of the abundance-dependent evolution of RSG spectral 
types to lower metallicities. We particularly note two 
RSGs with aberrantly late spectral types, which we con- 
sider to be robust variable candidates. Finally, we con- 
sider the implications of these new samples and the op- 
portunities they offer for future work (§ 4). 

2. DATA ACQUISITION 
2.1. Host Galaxies 

WLM was originally discovered by Wolf (1923) and 
Melotte (1926). It is classified as an Ir IV-V (van den 
Bergh 1966) and has a star formation rate (SFR) of 0.001 
Mq yr"i (Hodge & Miller 1995, Hunter & Elmegreen 
2004). Urbaneja et al. (2008) determine a distance mod- 
ulus to WLM of 24.99 ± 0.10 mag; Gieren et al. (2008) 
find a very similar distance of 24.92 ± 0.04 based on mul- 
tiwavelength Cepheid photometry. Massey et al. (2007b) 
obtained broad-band photometry for 7,656 stars seen to- 
wards WLM as part of the Local Group Galaxy Survey 
(LGGS), and determined a total E{B -V) = 0.07 ±0.05 
for this galaxy (foreground E{B — V) = 0.03 from 
Schlegel ct al. 1998). 

Lee et al. (2005) give a nebular abundance for WLM 
of log(0/H) + 12 = 7.83 ± 0.06. From an examina- 
tion of blue supergiants in WLM, Bresolin et al. (2006) 
find an identical abundance of log(0/H) -|- 12 = 7.83 it 
0.12. Venn et al. (2003) study two blue supergiants in 
WLM, and find that one of these stars has an extremely 
high oxygen abundance, with log(0/H) -|- 12 ~ 8.45; Ur- 
baneja et al. (2008) find similar results for this unusual 
star. Venn et al. (2003) speculate that this may be due 
to spatial chemical inhomogeneities in WLM; however. 



Lee et al. (2005) find no support for this based on their 
spectroscopic survey of HII regions in WLM. Finally, Ur- 
baneja ct al. (2008) also found comparable abundances 
to the Lee et al. (2005) nebular metallicity in their study 
of B and A supergiants in WLM. Venn et al. (2003), 
Lee et al. (2005), and Bresolin et al. (2006) all deter- 
mined specific elemental stellar abundances for Fe, Mg, 
N, and Si. Mg appears to be enhanced in WLM su- 
pergiants relative to what is expected from the nebular 
abundance; Bresolin et al. (2006) find enhanced N abun- 
dances in their blue supergiant spectra, but speculate 
that this is a consequence of enhanced rotational mix- 
ing effects, an effect which is predicted to be higher in 
low-metallicity environments (Macder & Meynet 2001, 
Leitherer 2008, Hirschi et al. 2008). Finally, Tramper 
et al. (2011) measure an unusually high mass loss rate 
for a WLM 0-star, and speculate that this unexpected 
result could potentially be attributed to unexpected ef- 
fects of low-metallicity evolution, such as wind clumping, 
pulsations, or magnetic fields. 

NGC 6822, originally discovered by Barnard (1885), 
is a barred dwarf irregular galaxy, also classified as an 
Ir IV-V (van den Bergh et al. 2000; Melotte 1926 notes 
the morphological similarities between NGC 6822 and 
WLM). This galaxy has an SFR of 0.01 Mq yr"! (Hunter 
& Elmegreen 2004). van den Bergh (2000) gives a dis- 
tance modulus to this galaxy of 23.45 ± 0.08 (0.50 Mpc); 
Pietrzyhski et al. (2004) find a slightly lower distance 
modulus of 23.34 ± 0.09. Massey et al. (2007b) obtained 
broad-band photometry for 51,877 stars seen towards 
NGC 6822 and calculate a total E{B -V) = 0.25 ± 0.02 
for this galaxy, a value which is dominated by the fore- 
ground E{B -V) = 0.22 (Schlegel ct al. 1998; it should 
also be noted than NGC 6822 is at a much lower galactic 
latitude than WLM, leading to a much higher number of 
foreground contaminants). 

Pagel et al. (1980) originally determined a nebular 
metallicity of log(0/H) + 12 = 8.25 ± 0.07 for NGC 
6822 based on observations of seven different HII regions. 
More recently, Lee et al. (2006) determine direct oxy- 
gen abundances from 5 HII regions and calculate a lower 
log(0/H) + 12 = 8.11 ± 0.1; however, this abundance 
is restricted to HII regions with detections of the auro- 
ral [OIII] A4363 line, suggesting that the survey may be 
biased towards lower-metallicity HII regions where the 
feature is strong enough to be detected (e.g. Garnett 
et al. 2004). Lee et al. (2006) also find no clear evi- 
dence for any abundance gradient in the galaxy. Abun- 
dance studies of individual supergiants in this host arc 
limited. Muschielok et al. (1999) find an overall [Fc/H] 
= -0.5 ± 0.2 dex based on spectra of three NGC 6822 B 
supergiants. Venn et al. (2001) determined a metallicity 
of log(0/H) + 12 = 8.37 ± 0.21 based on high-resolution 
spectroscopic studies of two A supergiants in NGC 6822, 
and find [Fe/H] = —0.49 ± 0.22, in agreement with the 
Muschielok et al. (1999) abundance and giving a metal- 
licity for NGC 6822 that is slightly higher than that of 
the SMC. However, most recent work in NGC 6822 has 
focused on the red giant and asymptotic giant branch 
stellar populations (e.g. Cioni & Habing 2005, Kang 
et al. 2006, Groenewegen et al. 2009); Cioni & Habing 
(2005) note that it is difficult to separate RSGs from the 
lower-mass giants in such work. 
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2.2. Sample Selection 

The RSGs in our samples were selected using photom- 
etry from LGGS. For both of these hosts, Massey et 
al. (2007b) originally plotted color-magnitude diagrams 
{B — V vs. V) and statistically "cleaned" the initial data 
to eliminate foreground dwarfs. Strong RSG populations 
were visible for both galaxies in this initial analysis. 

Massey (1998) originally identified a sample of RSGs 
in NGC 6822 (as well as samples in M31 and M33) us- 
ing the two-color method, applying a cut-off line in the 

V — R vs. B — V diagram that separates RSGs (above 
the line) from foreground dwarf contaminants (below the 
line) due to surface gravity effects. To select samples of 
RSG candidates from the LGGS data for our observa- 
tions, we began by selecting stars with V — R> 0.6 and 

V <20, restricting the sample to red stars bright enough 
to be observable with the Baade 6.5m telescope in expo- 
sure times of as little as ^2 hours. With these criteria 
in place, the initial sample selection was done using the 
Massey (1998) two-color method, applying a cutoff line 
of B -V ^ 1.25 X (y - R) + 0.45 in y - i? vs. B -V 
diagrams for both galaxies (Figure 1). For WLM this 
yielded an initial list of 18 potential RSG targets. Upon 
closer examination, one of these targets was found to be 
a galaxy in BVR images, and one more was removed due 
to crowding issues, yielding a final list of 16 RSGs to be 
observed in WLM. For NGC 6822, the V, V ~ R, and 
B — V criteria produced an initial list of 182 potential 
RSG targets. After removing 46 targets due to crowding 
issues, we were left with a final list of 136 RSG to observe 
in NGC 6822. 

Even with foreground dwarfs eliminated, red giants 
in the halo of the Milky Way can also potentially con- 
taminate samples of RSGs in Local Group galaxies (e.g. 
Levesque et al. 2007, Massey et al. 2009). For NGC 6822 
and WLM, we use the Besangon models of the Milky Way 
(Robin et al. 2003)0 to generate the expected number of 
stars as a function of magnitude and color in the square 
degree of sky centered on each galaxy. After applying the 
magnitude and color cuts used to select our RSG sample, 
we find only one remaining foreground contaminant star 
in each field, corresponding to a low likelihood of halo 
giant contamination in each case - 6% for the WLM and 
1% for NGC 6822. When applied to our final sample size 
in each galaxy, we therefore conclude that all of the stars 
selected here are bona fide RSG members of NGC 6822 
or WLM. 

2.3. Observations and Data Reduction 

The RSGs were observed using IMACS (Dressier et al. 
2011) on the Magellan Baade 6.5-meter telescope at Las 
Campanas Observatory on 07-08 August 2008. Condi- 
tions were clear on the first night, with sporadic high 
cirrus present on the second night; seeing on both nights 
was 0.5-0.6". The observations were obtained using five 
multi-object slitlet masks, three for WLM and two for 
NGC 6822, with slit widths of 1.2". These masks per- 
mitted us to observe all 16 WLM RSGs and 67 of the 
NGC 6822 RSGs. While all 16 WLM RSGs could have 
been included in two masks, three were advantageous as 
multiple orientations allowed us to account for the poten- 

® Available at |http: / /model. obs-besancon.fr/] 



tial effects of crowding. The IMACS f/4 camera gave us 
a field of view of 15' x 15', which is well-matched to the 
size of the galaxies. Using the 300 line mm~^ grating, we 
were able to obtain spectra with a resolution of lOA with 
an exposure time of 3 hours on each mask. The three 
WLM masks were observed at average airmasses of 1.05, 
1.25, and 1.38, while the two NGC 6822 masks were ob- 
served at average airmasses of 1.37 and 1.39. Due to the 
fixed orientation required for multislit masks, our obser- 
vations were not taken at the parallactic angle; however, 
IMACS has an Atmospheric Dispersion Corrector (ADC) 
that compensates for atmospheric differential refraction. 

We required continuous spectral coverage between 
^5000A and 8OOOA for our analyses. This initially posed 
a problem due to the configuration of the IMACS Mo- 
saicl detector at f/4: the detector is an 8192 x 8192 
CCD camera which consists of eight thinned 2K x 4K 
X 15^ CCD chips, arranged in a 2x4 grid with gaps 
of '^62 pixels (12.4") between the individual chips. With 
one spectrum extending across four of these chips, at any 
single grating tilt sections of the spectrum would fall into 
these gaps between the CCD chips (see Figure 2 for an 
example of our raw multislit spectroscopic data from Mo- 
saicl). This was avoided by observing each mask at two 
different grating tilts, 5.6 degrees (central wavelength of 
6OI5A) and 5.8 de grees (central wavelength of 6225 A), 
in order to ensure complete spectral coverage for each 
spectrum. 

We reduced the data using IRAF0, rather than the 
Carnegie Observatories System for MultiObject Spec- 
troscopy (COSMOS) data reduction package for IMACS. 
This permitted us to take a careful step-by-step approach 
to the delicate task of flux-calibrating multislit spectro- 
scopic observations across the individual CCD chips of 
the Mosaicl detector. We observed several spectropho- 
tometric standards for each mask (LTT1020, LTT7379, 
LTT9239, and CD-34 324; Hamuy et al. 1992, 1994), 
placing each standard in central slits of our masks on 
both the upper and lower rows of CCD chips and observ- 
ing each at both grating tilts. To ensure the best possible 
flux calibration, we then reduced the data for each indi- 
vidual chip separately, generating a specific sensitivity 
function (using IRAF's sensf unc task in the onedspec 
package) derived from the portion of the spcctrophoto- 
metric standard spectrum that fell on the chip. The RSG 
spectra on each chip were then calibrated using this sensi- 
tivity function and the calibrate function in IRAF, and 
finally the individual calibrated components of each RSG 
spectrum were median-combined using IRAF's scombine 
task to generate one complete stellar spectrum. We note 
that there are minor "plateaus" in the stellar spectra at 
-6000A-6300A and ~7600A-7800A, which correspond to 
the gaps in the CCD mosaic and do not impact our de- 
termination of the RSG spectral types. 

The main difficulty with flux calibrating the data in 
this manner arose when considering the horizontal (dis- 
persion axis) geometry of the multislit masks. With mul- 
tiple slits placed at a variety of positions relative to the 
grating, the wavelength coverage of each individual spec- 
trum varied slightly as a result. In the case of WLM, 

* IRAF is distributed by NOAO, which is operated by AURA, 
Inc., under cooperative agreement with the NSF. 



4 



Levesque & Massey 



where the geometry of the galaxy is narrow in the east- 
west dimension, this did not pose a significant problem 
- the wavelength coverage of the RSGs and the spec- 
trophotometric standards were comparable as a whole 
and allowed us to flux-calibrate the key 5000A-8000A re- 
gion in all of the spectra. However, for the more extended 
east-west geometry of NGC 6822 this effect was consider- 
able, and the spectral coverage of the spectrophotometric 
standards relative to the RSGs was insufficient. In Fig- 
ure 2 we present the raw data from Mosaic 1 for one of our 
NGC 6822 multislit masks; the bright sky lines illustrate 
the wide range in wavelength coverage spanned by the 
34 spectra observed with this mask. As a result, we were 
able to flux calibrate the WLM RSGs, but not the NGC 
6822 sample, which we leave in normalized relative flux 
units for the remainder of our analyses. For the WLM 
RSGs. 6 of the stars had poor S/N or poor wavelength 
coverage; however, we were able to successfully flux cal- 
ibrate the remaining II spectra. Based on our reduc- 
tion and analyses, we recommend that any future studies 
seeking to flux-calibrate multislit data should take par- 
ticular care to observe spectrophotometric standards at 
multiple positions along the dispersion axis of a slitmask. 

Once the flux calibration for our WLM RSGs was com- 
plete, we calculated approximate (monochromatic) V 
band magnitudes based on the stars' 5556A fluxes, and 
compared these to the V magnitudes determined from 
the LGGS data (Massey et al. 2007b). We found that 
all of our calibrations for each mask were consistently 
bright by uniform amounts, with standard deviations of 
0.1 mag. This suggested that the disparity we were see- 
ing was simply a grey shift in the spectrum. To con- 
flrm this we had to match the fluxes we determined for 
our flux-calibrated spectrophotometric standards, across 
a full range of bandwidths ranging from 5000A to 8OOOA. 
Comparing our measured fluxes to the standard star flies 
in IRAF's onedstds directory, we found that our flux- 
calibrated spectra agreed excellently with the measured 
values, differing by small uniform amounts with stan- 
dard deviations of several hundredths of a magnitude. 
Based on this, we can confidently claim that the dis- 
agreement between the V fluxes in our spectra and the 
V fluxes in the LGGS data are the result of an incorrect 
greyshift, rather than a problem with the flux calibration 
as a whole. We correct for this by manually greyshifting 
the observations to agree with the fluxes from LGGS. 
As a result, our observations cannot be used to deter- 
mine photometric data for these stars, and the spectral 
coverage and flux precision is insufficient for stellar at- 
mosphere model fitting. However, the overall shape of 
the stellar continuum has been properly determined and 
preserved, and can be used in comparisons with other 
RSG spectra to determine TiO band depths and spec- 
tral types. 

3. SPECTRAL TYPES 

Our reduced NGC 6822 and WLM RSG stellar spec- 
tra were compared to observations of spectral standards 
from Morgan & Kecnan (1973) by Levesque ct al. (2005). 
Since this work originally found that several of the spec- 
tral "standards" had to be reclassified, we also supple- 
ment our comparison spectra with additional RSGs from 
Levesque et al. (2005, 2006), to ensure consistency in our 



spectral types. 

For stars with late K or M spectral types, spectral 
types are based on the strengths of the TiO bands (which 
themselves are sensitive to physical properties such as 
temperature and metallicity) , with stronger bands corre- 
sponding to later spectral types. Specifically, the types 
are primarily determined from the 6I58A, 6658A, and 
7054A bands following Jaschek & Jaschek (1990), with 
TiO bands further in the blue (5167A, 5448A, 5847A) 
serving as secondary confirmations of the quality of the 
flt. 

For early-K type stars, this classification is consider- 
ably more challenging. Typically, we base these spectral 
types on the strength of the G band and the Ca I A4226 
line, since at this wavelength coverage and spectral res- 
olution (~I0A) the spectrum of a K-type star is nearly 
featureless. However, these features are not included in 
our more limited spectral coverage for these RSGs, and 
we must therefore base our early-K spectral types en- 
tirely on the overall goodness of the spectral continuum 
fit. 

For the NGC 6822 RSGs, we normalized our compari- 
son star spectra for use with the non-flux-calibrated data. 
This allowed us to determine a best flt for stars with 
good S/N and evidence of TiO bands in their spectra; 
however, this also imposes a limitation on determining 
conclusive spectral types for NGC 6822 RSGs with early- 
and mid-K spectral types, as we cannot base our flts on 
the overall SED shape. We must also exclude a number 
of RSGs with poor or noisy spectra where emerging TiO 
features cannot be confldently distinguished. As a result, 
we have limited our NGC 6822 spectral types to stars 
where we can confidently assign a spectral type of K5 I 
or later, the earliest type where TiO features begin to 
emerge. This results in a total sample of 16 NGC 6822 
RSGs, with a strong bias towards the late-type RSGs. 
In the case of the WLM RSGs, we were able to flt all 
II of our flux-calibrated spectra with the spectral stan- 
dards, and could extend our spectral type determinations 
to the early-K subtypes; however, it should be stressed 
that these K-type classifications have a lower confidence 
that the TiO-based late-K and M spectral types. While 
we can confidently distinguish between "early" K type 
stars with no visible TiO features (KO-I I or K2-3 I) and 
late-K types where TiO features begin to emerge (K5 I), 
the distinction between these early-type classifications is 
subject to the overall effects of reddening and surface 
gravity on the spectrum, as well as the goodness of the 
flux calibration. 

The spectral types and colors for these stars are given 
in Table I ; color indices were determined using B VR pho- 
tometry from LGGS and Kg photometry from 2MASS, 
while previous RSG identiflcations and spectral types are 
taken from Massey (1998) or targets from Humphreys 
(1980) cross-referenced by Massey et al. (2007b). It 
should be noted that, while a relation has previously 
been seen between {V — K)o and Toff (Levesque et al. 
2005, 2006), we cannot use this relation with our cur- 
rent data. Converting from V — Kg to {V — K)o, where 
{V - K)o- = V-iKs + 0.04) - (0.88 x Ay) following 
Schlegel et al. (1998) and Carpenter (2001), requires a de- 
termination of Av for each star. At present we can only 
give conservative lower limits of Ay = 0.09 for WLM 
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and Ay = 0.68 for NGC 6822 based on foreground ex- 
tinction (Massey ot al. 2007b). It is well-known that 
RSGs can display substantial amounts of excess redden- 
ing due to circumstellar dust (e.g. Massey et al. 2005). 
In addition, both Lee et al. (2005) and Urbaneja et al. 
(2008) have noted evidence of variations in reddening 
throughout WLM in particular; Urbaneja ct al. (2008) 
find Ay values ranging from 0.09 to 0.90 in their B 
and A supergiant sample. With the uncertain redden- 
ing for each RSG in our sample, we are therefore not 
currently equipped to estimate Tes for the sample. Simi- 
larly, determining bolometric luminosities is also not pos- 
sible with the ciirrcnt data since, in addition to account- 
ing for imknown reddening effects, RSGs have significant 
Teff-dependent bolometric corrections (a 10% error in T^g 
corresponds to a factor of 2 error in bolometric luminos- 
ity; Kurucz 1992, Massey & Olsen 2003). In addition to 
these RSGs, in Table 2 we Ust the additional NGC 6822 
RSGs from our sample whose spectral types could not be 
determined due to weak TiO features or noisy spectra. 

The spectra of RSGs that we were able to assign spec- 
tral types to are shown in Figures 3 (WLM) and 4 (NGC 
6822) with spectral coverage from ~5000A - 8OOOA, From 
these results, two RSGs immediately stand out in our 
samples as having unusually late spectral types relative 
to the overall RSG population in their host galaxies: 
J000158.14-152332.2 in WLM, with a spectral type of 
M 3 I, and J194453.46-144552.6 in NGC 6822, with a 
spectral type of M4.5 I. Both of these stars have spectra 
dominated by extremely strong TiO features. Interest- 
ingly, the spectrum of J000158. 14-152332. 2 also shows 
Ha in emission, and the SED shows evidence of substan- 
tial reddening compared to the other WLM RSGs. Both 
of these properties, along with its relatively late spectral 
subtype, are also similar to the unusual late- type dust- 
enshrouded RSG WOH G64 in the LMC (Ohnaka et al. 
2008, Levesque et al. 2009). From our test of halo giant 
contamination using the Besancon models and our elim- 
ination of foreground dwarfs using the Massey (1998) 
two-color method, we feel confident that both of these 
are extragalactic RSGs, 

The extreme nature of these stars' spectral types, and 
their significance in studies of unusual RSGs, is illus- 
trated in Figure 5. Here we plot histograms of spectral 
types for RSGs from the Milky Way (Levesque et al. 
2005), the Large and Small Magellanic Clouds (Levesque 
et al. 2006, 2007, 2009; Massey et al. 2007b), and NGC 
6822 and WLM (this work). It should be noted that the 
selection criteria for RSGs across different host galaxies 
results in a variation of the parameter space that is sam- 
pled; while the Milky Way RSGs were selected based on 
membership in OB associations, the extragalactic RSG 
samples were selected based on their V magnitudes and 
colors. In addition, more observations of RSGs in WLM 
and NGC 6822 are required to improve the complete- 
ness of the plot (particularly for early-K supergiants in 
NGC 6822). Despite these deficiencies, we can clearly see 
that these samples follow the trend originally described 
by EUas et al. (1985) and Massey & Olsen (2003), with 
the average spectral subtype for RSGs shifting to ear- 
lier types at lower metallicities as a result of metallicity 
effects on TiO band strengths and the effective temper- 
ature of the Hayashi limit. 



In spite of this overall trend, and the expected accom- 
panying restrictions of the Hayashi limit, unusually late- 
M supergiants are still present in all of the low-metallicity 
samples, J000158,14-152332.2 and J194453.46-144552.6 
are both marked as clear outliers when compared to 
the other RSGs in their host galaxies. Also labeled are 
the variable late-type outliers in the Magellanic Clouds 
(LMC 170452, SMC 046662, SMC 055188, and SMC 
050028; Massey et al. 2007, Levesque et al. 2007) and the 
unusual dust-enshrouded RSG WOH G64 in the LMC. 
From this histogram, it is clear that both J000158.14- 
152332.2 and J194453.46-144552.6 are excellent candi- 
dates for follow-up studies searching for variable or diist- 
enshrouded RSGs in these low-metallicity host galaxies. 

In addition, there is one RSG in our NGC 6822 sample, 
J194449. 96-144333.5, whose spectral type was found to 
be at odds with previous determinations. Massey (1998) 
assigned this star a spectral type of M2.5, while we as- 
sign it a spectral type of K5. In Figure 6 we show the 
spectnun used in the Massey (1998) classification, taken 
on 8 September 1997 at the Kitt Peak National Obser- 
vatory 4-meter telescope, and compare it to the spec- 
trum from this work. From these observations, it is clear 
that the star's spectral type appears to be variable, with 
the 1997 spectrum showing much stronger TiO bands. 
The 1997 spectrum also shows an unusual strong Ha 
absorption line (equivalent width ~6A) that is weaker 
in the 2008 spectrum (equivalent width ~lA). For com- 
parison, Figure 7 shows three different spectra of the 
previously-identified variable HV 11423, an RSG in the 
SMC. These spectra show very similar variations to our 
observations of J194449. 96-144333. 5; the star appears as 
a KO I in 2004, an M4 I in 2005, and an M2.5 I in new 
observations of the star taken with IMACS on Magellan 
in August of 2010. We therefore suggest that, in ad- 
dition to J000158. 14-152332. 2 and J194453.46-144552.6, 
J194449, 96-144333. 5 is a prime candidate for follow-up 
studies examining variable behavior. In addition, this 
star illustrates that type- variable RSGs cannot necessar- 
ily be identified by aberrantly late spectral types in a 
single epoch of observations, since it is possible that they 
were simply captured in a early-type state. 

4. DISCUSSION AND FUTURE WORK 

We observed extragalactic RSG populations in two 
low-metallicity Local Group galaxies using IMACS. By 
combining these multislit observations with careful re- 
ductions and flux calibrations, we have acquired spectro- 
scopic data for 16 RSGs in NGC 6822 and spectropho- 
tometric data for 11 RSGs in WLM. We have deter- 
mined spectral types for each of our RSG spectra, and 
use these data to extend the metallicity-dependent evo- 
lution of average RSG spectral types with abundance 
to new and lower metallicities. In addition, we have un- 
covered two RSGs, .1194453,46-144552,6 and J000158.14- 
152332,2, with anomalously late spectral types relative 
to their host populations, and a third RSG, J194449.96- 
144333.5, that shows evidence of variable behavior from 
past spectra. 

Understanding the impact of metallicity on the ob- 
served properties of RSGs is vital to studies of extra- 
galactic massive stars. The physical properties of RSGs, 
along with the phenomena driving their mass loss pro- 
cesses, are critical in local studies of supernova progen- 
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itors (e.g. Smartt ct al. 2009, Gcorgy 2012) and in 
high-redshift studies of yoking stellar populations where 
RSGs are the dominant dust producers (Massey et al. 
2005). In addition, the trend toward earlier spectral 
types poses a potential problem for fully sampling the 
evolved supergiant populations in lower-metallicity hosts 
such as NGC 6822, WLM, and the SMC. Any such survey 
must make an effort to include yellow supergiants, due 
to the metallicity-dependent evolutionary effects on the 
Hayashi limit and the resulting skew towards an evolved 
supergiant population with slightly bluer colors at lower 
metallicities as a whole. However, removing foreground 
contaminants in yellow supergiant samples is challenging 
and requires detailed kinematic treatments (e.g. Drout 
et al. 2009, Neugent et al. 2010). 

Finally, the potential addition of two new stars to the 
sample of variable RSGs is extremely valuable. This is 
a small and poorly-understood sample of evolved mas- 
sive stars, one that would benefit from additional obser- 
vations or regular monitoring programs. New candidates 
could be detected through photometric observations, due 
to variations in these stars' color indices, but unambigu- 
ous detections and confirmation of variability in stars 
such as J194449.96-144333.5, J194453.46-144552.6, and 
J000158. 14-152332. 2 will require multiple spectroscopic 



observations. Confirming variable behavior in new can- 
didates, and beginning to understand the physical phe- 
nomena driving these stars' behavior, will require addi- 
tional observations spanning a baseline of several years. 

We thank the anonymous referee for their valuable 
feedback on this manuscript. We also gratefully ac- 
knowledge valuable discussions with Nidia Morrell. We 
are thankful to the staff of Las Campanas Observatory 
for their assistance and hospitality during our observa- 
tions. Magellan telescope time was granted by NO AO, 
through the Telescope System Instrumentation Program 
(TSIP). TSIP is funded by NSF. This paper made use 
of data from the Two Micron AU Sky Survey (2MASS), 
which is a joint project of the University of Massachusetts 
and the Infrared Processes and Analysis Center, Cali- 
fornia Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National 
Science Foundation. E. M. L. is supported by NASA 
through Einstein Postdoctoral Fellowship grant num- 
ber PFO-110075 awarded by the Chandra X-ray Center, 
which is operated by the Smithsonian Astrophysical Ob- 
servatory for NASA under contract NAS8-03060. This 
work was supported by the National Science Foundation 
through AST-1008020 to P. M. 



REFERENCES 



Dressier, A., et al. 2011, PASP, 123, 288 
Barnard, E. E. 1884, AN, 110, 125 

Blair, W. P., Kirshner, R. P., & Chevalier, R. A. 1982, ApJ, 254, 
50 

Bresolin, P., Pietrzynski, G., Urbaneja, M. A., Gieren, W., 

Kudritzki, R.-P, & Venn, K. A. 2006, ApJ, 648, 1007 
Carpenter, J. M. 2001, AJ, 121, 2851 
Cioni, M.-R. L. & Habing, H. L. 2005, A&A, 429, 837 
Charbonnel, C, Mcynet, C, Maeder, A., Schaller, G., & Schaerer, 

D. 1993, A&AS, 101, 415 
Drout, M. R., Massey, P., Meynet, G., Tokarz, S., & Caldwell, N. 

2009, ApJ, 703, 441 
Elias, J. H., Frogel, J. A., & Humphreys, R. M. 1985, ApJS, 57, 91 
Garnett, D. R., Kennicutt, R. C, & Bresohn, F. 2004, ApJ, 607, 

L21 

Georgy, C. 2012, A&A, 538, L8 
Gieren, W., et al. 2008, ApJ, 683, 611 

Groenewegen, M. A. T., Langon, A., & Marescaux, M. 2009, A&A, 
504, 1031 

Gustafsson, B., Edvardsson, B., Eriksson, K., Mizumo-Wiedner, 
M., Jorgensen, U. G., & Plcz, B. 2003, in ASP Conf. Scr. 288, 
Stellar Atmosphere Modeling, ed. I Hubcny, D. Mihalas, & K. 
Werner (San Francisco; ASP), 331 

Gustafsson, B., Edvardsson, B., Eriksson, K., Jorgensen, U. G., 
Nordlund, A, & Plez, B. 2008, A&A, 486, 951 

Hamuy, M., Suntzeff, N. B., Heathcotc, S. R., Walker, A. R., 
Gigoux, P., & Phillips, M. M. 1994, PASP, 106, 566 

Hamuy, M., Walker, A. R., Suntzeff, N. B., Gigoux, P., Heathcote, 
S. R., Phillips, M. M. 1992, PASP, 104, 533 

Hayashi, C. & Hoshi, R. 1961, PASJ, 13, 552 

Hirsclii, R., Chiappini, C., Meynet, G., Maeder, A., & Ekstrom, S. 

2008, in lAU Symp. 250, Massive Stars as Cosmic Engines, ed. 

F. Bresolin, P. A. Crowther, & J. Puis (Cambridge: Cambridge 

Univ. Press), 217 
Hodge, P. W. & Miller, B. W. 1995, ApJ, 451, 176 
Hunter, D. A. & Elmegreen, B. G. 2004, AJ, 128, 2170 
Jaschek, C. & Jaschek, M. 1990, The Classification of Stars 

(Cambridge: Cambridge Univ. Press) 
Kang, A., et al. 2006, A&A, 454, 717 

Kurucz, R. L. 1992, in lAU Symp. 149, The Stellar Populations of 
Galaxies, ed. B. Barbuy & A. Renzini (Dordrecht: Kluwer), 225 

Lee, H., Skillman, E. D., Cannon, J. M., Ja<;kson, D. C, Gehrz, R. 
D., Polomski, E. F., & Woodward, C. E. 2006, ApJ, 647, 970 

Lee, H., Skillman, E. D., & Venn, K. A. 2005, ApJ, 620, 223 



Leitherer, C. 2008, in lAU Symp 255, Low-Metallicity Star 
Formation: From the First Stars to Dwarf Galaxies, ed. L. K. 
Hunt, S. Madden, & R. Schneider (Cambridge: CUP), 305 

Levesque, E. M., Massey, P., Olsen, K. A. G., & Plez. B. 2007, 
ApJ, 667, 202 

Levesque, E. M., Massey, P., Olsen, K.A.G., Plez, B., Josselin, E., 

Maeder, A., Meynet, G 2005, ApJ, 628, 973 
Levesque, E. M., Massey, P., Olsen, K. A. G., Plez, B., Meynet, 

G., & Maeder, A. 2006, ApJ, 645, 1102 
Maeder, A. & Meynet, G. 2001, A&A, 373, 555 
Massey P. 1998, ApJ, 501, 153 
Massey P. 2003, ARA&A, 41, 15 

Massey, P., Levesque, E. M., Olsen, K. A. G., Plez, B., & Skiff, B. 

A. 2007a, 660, 301 
Massey P. & Olsen, K. A. G. 2003, AJ, 126, 2867 
Massey, P., Olsen, K. A. G., Hodge, P. W., Jacoby G. H., McNeill, 

R. T., Smith, R. C, & Strong, S. B. 2007b, AJ, 133, 2393 
Massey, P., Plez, B., Levesque, E. M., & Olsen, K. A. G. 2005, 

ApJ, 634, 1286 

Massey, P., Silva, D. R., Levesque, E. M., Plez, B., & Clayton, G. 

C. 2009, ApJ, 703, 420 
Melotte, P. J. 1926, MNRAS, 86, 636 
Meynet, G. & Maeder, A. 2005, A&A, 429, 581 
Morgan, W. W. & Keenan, P. C. 1973, ARA&A, 365, 347 
Muschielok, B., et al. 1999, A&A, 352, 40 

Neugent, K. F., Massey P., Skiff, B., Drout, M., Meynet, G., & 

Olsen, K. A. G. 2010, ApJ, 719, 1784 
Ohnaka, K., Driebe, T., Hofmann, K.-H., Weigelt, G., & 

WIttkowski, M. 2008, A&A, 484, 371 
Pagel, B. E. J., Edmunds, M. G., & Smith, G. 1980, MNRAS, 193, 

219 

Pietrzynski, G., Gieren, W., Udalski, A., Bresolin, F., Kudritzki, 
R. P., Soszyhski, I., Szymahski, M., & Kubiak, M. 2004, AJ, 
128, 2815 

Plez, B. 2003, in ASP Conf. Ser. 298, GAIA Spectroscopy; Science 
and Technology, ed. U. Munari (San Francisco: ASP), 189 

Robin, A. C, Reyle, Ci, Derricre, S., & Picaud, S. 2003, A&A, 
409, 523 

Schaerer, D., Meynet, G., Maeder, A., & Schaller, G. 1993, A&AS, 
98, 523 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Smartt, S. J. Eldridge, J. J., Crockett, R. M., & Maund, J. R. 2009, 

MNRAS, 395, 1409 
Tramper, F., Sana, H., do Koter, A., & Kaper, L. 2011, ApJ, 741, 

L8 



RSGs in NGC 6822 and WLM 



7 



TABLE 1 

WLM AND NGC 6822 Red Supergiants 



Star (LGGS) 


Type 




B 


V - R'' 


V - Ks"^'^^ 


Previous Work 












i'ype 


Hef^ 


WLM 


J000153. 17-152813.4 


KO-1 


19.34 


1.66 


0.83 


4.07 






J000156. 77-152839. 6 


K2-3 


17.61 


1.94 


0.99 


3.81 






J000156. 87-153122.3 


KO-1 


18.76 


1.64 


0.85 


3.19 






J000157.01-152954.0 


KO-1 


18.70 


1.64 


0.85 


3.43 






J000157.55-152915.8 


KO-1 


19.27 


1.64 


0.86 


3.55 






J000157.96-152803.1 


KO-1 


19.80 


1.53 


0.77 








J000158.14-152332.2 


M3 


19.61 


1.72 


0.99 








J000158.74-152245.5 


KO-1 


18.46 


1.68 


0.90 


3.12 






J000159.61-153059.9 


K2-3 


18.98 


1.78 


0.91 


3.69 






J000200.81-153115.7 


KO-1 


18.69 


1.78 


0.91 


3.95 






J000203.04-153033.7 


K5 


18.68 


1.98 


1.04 


4.41 






NGC 6822 


J194445.76-145221.2 


Ml 


17.71 


2.22 


1.21 


4.87 






J194447.56-145215.4 


K5 


19.13 


2.05 


1.12 


4.80 






J194447.81-145052.5 


Ml 


18.51 


2.23 


1.29 


5.26 






J194449. 96-144333. 5 


K5 


18.07 


2.21 


1.29 


5.24 


M2.5 I 




J194450. 12-144637.9 


K5 


19.84 


1.71 


0.90 








J194450. 44-144410.0 


M2 


18.46 


2.00 


1.05 


4.29 






J194453. 46-144540.1 


K5 


19.14 


1.74 


0.97 


4.27 






J194453. 46-144552. 6 


M4.5 


18.43 


1.93 


1.05 


4.35 






J194453.96-144424.3 


K5 


18.19 


2.01 


1.04 


4.46 




1 


J194454.46-144806.2 


Ml 


18.56 


2.00 


1.21 


5.35 






J194454.54-145127.1 


MO 


17.05 


2.25 


1.19 


4.71 


MO-1 I 


2 


J194455.70-145155.4 


MO 


16.91 


2.20 


1.17 


4.55 


Ml-2 I 


2 


J194455.93-144719.6 


K5 


19.56 


1.87 


1.04 


4.46 






J194457.31-144920.2 


Ml 


17.41 


2.28 


1.21 


4.93 


Ml I 


2 


J194459.86-144515.4 


Ml 


16.93 


2.00 


1.00 


4.27 


MO I 


1,2 


J194503.58-144337.6 


MO 


19.23 


1.91 


1.01 


4.27 







^ From Massoy ot al. (2007b). 

^ Ks magnitudes are taken from 2MASS wiiere available. 

References: 1) Massey (1998). 2) Spectral type from Humphreys (1980), cross-identifications 
from Massey et al. (2007b). 
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(Cambridge: Cambridge Univ. Press) 
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TABLE 2 

Additional NGC 6822 Red Supergiants 



Star (LGGS) 






B 


- 


V 


- 


V - 




J 194449 


62-145327 


3 


19 


470 


1 


769 





957 


4 


184 


J 194451 


7o-144d1d 


7 


19 


848 


1 


720 





895 






J 194400 


1 o -\ A tz A r\r- 

1o-1454Ud 


5 


19 


604 


1 


839 





982 


4 


375 


J 194453 


66-144617 


9 


19 


050 


2 


030 


1 


169 


4 


872 


J 194450 


Do-145z45 


8 


19 


252 


1 


816 





979 


4 


270 


J 194454 


01 T A A A '~7 

5l-144o47 


8 


18 


050 


2 


239 


1 


284 


5 


324 


T 1 A A A 

J 194454 


5d-145z1d 


3 


18 


768 


1 


854 


1 


062 






J 194455 


02-144607 


1 


19 


465 


1 


900 





993 


4 


200 


J 194455 


'\ (\ ^ A A ^ A A 

19-144344 


8 


19 


814 


1 


652 





931 


4 


845 


J 194455 


37-145229 


2 


19 


272 


1 


798 





976 


3 


941 


T 1 f\ A AV ^ 

J 194455 


97-145109 


4 


19 


401 


1 


840 


1 


013 


4 


302 


J 194450 


77-144927 





19 


879 


1 


655 





883 






J 194457 


10-145141 


1 


19 


950 


1 


700 





961 






J 194457 


38-145408 





18 


508 


2 


095 


1 


114 


4 


555 


T 1 Ct A A fZ'y 

J 194457 


69-144953 


1 


19 


137 


1 


857 


1 


005 


4 


594 


J 194457 


nr* 1 A A o A n' 

9d- 144847 


9 


19 


991 


1 


760 





946 




_■ 


J 194450 


ir* -\ A A AC r- 

1d-14445d 





18 


329 


2 


067 


1 


126 


4 


751 


J 194458 


17-144527 


8 


18 


143 


1 


924 





977 






J 194458 


45-14545o 


8 


19 


553 


1 


580 





836 






J 194458 


64-1450U7 


4 


19 


669 


1 


733 





912 






J 194459 


OO 1 A ACltZ fZ 

22-144955 





19 


764 


1 


730 





940 






J 194459 


A'X 1 /I cri f\f 

41-1451UO 


8 


19 


877 


1 


721 





948 






J 194500 


1 r\ ^ A A r\r\ "1 

10-144901 


7 


19 


789 


1 


664 





930 


4 


346 


J 194500 


1 A A ^ f\ A 

32-144104 


8 


18 


246 


1 


999 


1 


062 


4 


553 


J 194500 


40-145333 


3 


19 


674 


1 


896 





998 


4 


158 


J 194500 


53-144820 


5 


18 


750 


2 


092 


1 


238 


5 


216 


J 194500 


rti Id coo A 

yi-145234 


3 


19 


450 


1 


788 





976 


4 


090 


J194501 


12-145437 


7 


18 


979 


1 


813 





971 


4 


239 


J194503 


27-144224 


5 


19 


901 


1 


718 





927 






J194503 


87-144247 


6 


19 


957 


2 


060 


1 


120 


4 


801 


J194504 


29-145406 


2 


19 


456 


1 


915 


1 


059 


4 


571 


J194504 


74-145422 


8 


19 


862 


1 


706 





926 


4 


574 


J194505 


74-145547 


6 


18 


862 


1 


892 


1 


021 


4 


361 


J194505 


84-145517 


3 


18 


935 


1 


774 





967 


4 


281 


J194506 


34-145537 





19 


382 


1 


951 


1 


041 


4 


283 


J194506 


35-145135 


5 


19 


155 


2 


076 


1 


140 


4 


672 


J 194507 


48-145107 


1 


19 


677 


1 


881 


1 


002 






J194507 


87-145435 


9 


19 


822 


1 


813 





943 






J194509 


55-145344 


6 


18 


838 


1 


880 





985 


4 


189 


J194515 


10-144803 


9 


19 


263 


1 


761 





781 


3 


919 



^ From Massey et al. (2007b). 

^ Kg magnitudes are taken from 2MASS where available. 
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Fig. 1. — V — R vs. B — V plots used to eliminate foreground dwarfs from LGGS photometry centered on WLM (left) and NGC 6822 
(right). Presumed RSGs with V < 20 are shown as red crosses; foreground dwarfs are plotted as black crosses. Targets with observations 
presented in this paper are overplotted as blue circles. 
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Fig. 2. — Reconstructed example of the IMACS Mosaicl CCD chip geometry, using raw data from our observatiori@0&f Sane of our two 
multislit masks for NGC 6822. Wavelength decreases from left to right in this image. The gaps between the individual CCD chips are 
shown, illustrating the need for observations at two different grating tilts to ensure complete spectral coverage. The varying x-axis positions 
of the bright sky lines illustrate the wide range in wavelength coverage for the 34 RSG spectra included on these masks. The chips in the 
bluest quarter of the CCD (far right) were found to have extremely poor S/N and were not used in our subsequent reductions and analyses. 
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Fig. 3. — Spectrophotometry and spectral types for our 11 WLM RSGs. The strong feature at 7600 A is the telluric A band. 
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5000 5500 
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Fig. 4. — Normalized spectroscopy and spectral types for our 16 NGC 6822 RSGs. The strong feature at 7600 A is the telluric A band. 
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Flc;. 5. - Histograms of RSG spectral types found in five different Local Group galaxies, plotted from top to bottom in order of decreasing 
metallicity. Data for the Milky Way RSGs arc taken from Levesque ct al. (2005); for the Magellanic Clouds spectral type data comes 
from Levesque ct al. (2006, 2007) and Massey et al. (2007b). Spectral types for NGC 6822 and WLM RSGs are from this work. In this 
comparison we can observe the progression of the dominant spectral type towards earlier types at lower metallicities that was previously 
discussed by Elias et al. (1985) and Massey & Olsen (2003). Late-type RSGs from Levesque et al. (2007) and this work are labeled. 
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Fig. 6.— Normalized spectra of the NGC 6822 RSG J194449.96-144333.5 from 1997 observations described in Massey (1998; top) and 
this work (bottom). The 1997 spectrum shows strong 5448A, 5847A, and 6158A TiO absorption bands, consistent with a spectral type of 
M2.5, while the 2008 spectrum shows only the weak 6158A TiO band, consistent with a K5 spectral type. 
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Fig. 7. — Moderate-resolution spectra of the unusual variable SMC 050028 in 2004 (top), 2005 (center), and 2010 (bottom). Observed 
spectra are shown in black; MARCS stellar atmosphere models (Plez 2003; Gustafsson et al. 2003, 2008) at SMC metallicity are shown in 
red. The variation in the star's spectral type and T^ff is apparent from the changing strengths of the TiO bands; using the MARCS models 
we determine a Toff of 4300 K from the 2004 data, 3500 K in 2005, and 3625 K in 2010. 



